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Guidance notes 


This computer book contains those sections of the chapters in Block A 
which require you to use Mathcad. Each of these chapters contains 
instructions as to when you should first refer to particular material in this 
computer book, so you are advised not to work on the activities here until 
you have reached the appropriate points in the chapters. 


In order to use this computer book, you will need the following Mathcad 


files. 


Chapter Al 


221A1-01 Fibonacci numbers 
221A1-02 Linear second-order recurrence sequences 
221A1-03 A Fibonacci sunflower (Optional) 


Chapter A2 


221A2-01 Ellipses in parametric form 
221A2-02 Hyperbolas in parametric form 
221A2-03 Focus, directrix and eccentricity (Optional) 


Chapter A3 
221A3-01 Isometries and triangles 


221A3-02 Isometries and conics 
221A3-03 Surface and contour plots (Optional) 


Instructions for installing these files onto your computer’s hard disk, and 
for opening them, are given in Chapter AO of MST121. 


The computer activities for each chapter also require you to work with 
Mathcad documents which you have created yourself. 


Activities based on software vary both in nature and in length. Sometimes 
the instructions for an activity appear only in the computer book; in other 
cases, instructions are given in the computer book and on screen. 
Feedback on an activity is sometimes provided on screen and sometimes 
given in the computer book. 


For advice on how each computer session fits into suggested study 
patterns, refer to the Study guides in the chapters. 


Chapter Al, Section 4 


Exploring linear second-order recurrence 
sequences with the computer 


Mathcad can quickly produce tables of terms of linear second-order 
recurrence sequences, and of expressions related to such sequences. These 
can help you to spot patterns in the sequences, and so form conjectures 
about them. In this section you will be introduced to the Mathcad 
techniques needed to produce such tables. You will then use examples of 
these tables to look for patterns and to form conjectures. 


The section ends with an optional subsection in which you are invited to 
explore the relationship between the pattern of florets on a sunflower head 
and the Fibonacci numbers, using a Mathcad file that plots a sunflower 
pattern. 


Further information on Mathcad features can be found in the MS221 
reference manual A Guide to Mathcad. 


4.1 Computer exploration of Fibonacci numbers 


In the first activity, you will see how to use Mathcad to calculate terms of 
the Fibonacci sequence F), using the recurrence system 


eee + FF, (n= 0,1,2,...), (4.1) See Chapter Al, Section 2. 


and how to display these terms in a convenient table. 


Activity 4.1 The Fibonacci recurrence system 


Open Mathcad file 221A1-01 Fibonacci numbers. Page 2 of the 
worksheet describes some key features of MS221 Mathcad worksheets. The Mathcad worksheets for 


M$221 have a similar design 
When you have finished with page 2, move to page 3 of the worksheet and to those for MST121. 


carry out Task 1. Remember to make your own 


A solution is given on page 4 of the worksheet. working copies of them. 


Comment 


© Range variables and subscripted variables play essential roles in 
calculating terms of a recurrence sequence in Mathcad. At the 
beginning of Task 1, you defined the range variable n as follows: 


n:=0,1..N —2. 


This ensures that the subsequent definitions of the subscripted 
variables Fo, F, and Fi,+9, 


Fo _— 0, F, _— : Fuse _— : ae = : 


together define all the Fibonacci numbers Fo, Fi,..., ny, and no 
others. Mathcad carries out the final definition here N — 1 times, once 


for each of the values 0,1,...,.N — 2 in the range of n. As it does so, : : 5 : : = : 
the subscript n + 2 takes each of the values in the range 2,3,...,N in 
turn. 


If you increase the value of N, 
remember to set it to 20 
again before continuing. 


See Chapter Al, Section 3. 
Recall that 


¢ = $(1+ v5), 
v= i(1— v5). 


You should still be working 
with Mathcad file 221A1-01. 
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© You displayed the Fibonacci numbers Fo, F,,..., Fy by entering ‘fF’ =’. 
This causes all the values of the subscripted variable to appear in a 
table. The table displays the subscripts on the left (if there are more 
than nine values), and scrolls if there are many values to display. 


You defined the variable N to be 20, to calculate the Fibonacci 
numbers Fo, F,,..., Fo9. You could calculate more Fibonacci numbers, 
by increasing the value of N. Mathcad automatically updates the 
values of all other variables whose definitions appear subsequently and 
depend on JN, for example, the range variable n := 0,1.. N — 2. 


© 


In Activity 4.1, you saw how to calculate and display terms of a sequence 
defined by a recurrence system. In the next activity, you will see a simple 
way to calculate and display terms of a sequence defined by a formula 
involving a variable such as n. This method will be used to produce tables 
showing the sequences defined by Binet’s formula, 


P= Ae ae oS 


and by Binet’s approximation, $”/V/5. 


Activity 4.2 Binet’s formula and Binet’s approximation 


Move to page 5 of the Mathcad worksheet, and carry out Task 2. Compare 
the values calculated using Binet’s formula and Binet’s approximation with 
the Fibonacci numbers Fo, F,, F2,...,F9, which were calculated using the 
recurrence system (4.1). 


Comment 


© The table for Binet’s formula shows that the formula does indeed give 
the same value for F,, as the recurrence system, for n = 0,1,..., 20. 


The table for Binet’s approximation makes it easy to see that F;, is 
indeed the nearest integer to the value given by Binet’s approximation, 
for n = 0,1,...,20. You can also see that, as n increases through the 
range 0,1,...,20, the values of Binet’s approximation alternate above 
and below the Fibonacci numbers, and become progressively closer to 
them. 


The tables of values on page 5 of the Mathcad worksheet are produced 
by evaluating formulas involving a range variable which has been 
defined previously. Each table displays the values taken by the 
corresponding formula as the range variable takes each value in its 
range. By default, the table scrolls if there are many values to display, 
but it can be resized to see all of the values without scrolling. 


© 


Here, the range variable is n; its definition, n := 0,1..N, appears after 
the words ‘Display range’. The first table is for the formula n and 
therefore displays 0,1,...,N. The second is for F;,, so this table 
displays Fo, F,,..., fn, which were all defined on page 3 of the 
worksheet. The other two tables are for Binet’s formula and Binet’s 
approximation. 


SECTION 4 Exploring linear second-order recurrence sequences with the computer 


Mathcad notes 


© Mathcad variables can be defined more than once, so they can take 
different values at different places in a worksheet. For example, the 
range variable n is defined as n := 0,1.. N — 2 on page 3 (to calculate 
the Fibonacci numbers) and as n := 0,1...N on page 5 (to display 
them). 


Each definition of a variable is used by Mathcad for all expressions 
involving the variable which appear below that definition in the 
worksheet, up to the point where a new definition of the same variable 
appears or the worksheet ends. 


© The table of Fibonacci numbers on page 5 of the worksheet appears 
different in style from that on page 4. The table obtained by entering 
‘F,, =’ on page 5 is displayed below the expression, whereas that 
obtained by entering ‘F' =’ on page 4 is displayed to the right, with an 
additional grey column to denote the subscript values. These 
differences (between tables containing the same output values) are due 
to different Mathcad display defaults being triggered in either case. 
The table characteristics can be altered by clicking on a value in the 
table with the right mouse button, and then choosing either 
‘Properties...’ or ‘Alignment’ from the resulting mini-menu. 


© ‘To enter a Greek letter in Mathcad, you can either click on the 
appropriate button on the ‘Greek’ toolbar, or type the equivalent 
Roman letter and then press [Ctrl]g. For example, if you type a 
followed by [Ctrl]g, then the ‘a’ will change into an ‘a’ (alpha). 


Now close Mathcad file 221A1-01. 


4.2 Exploring patterns in linear second-order 
recurrence sequences 


In the next series of activities, you are invited to explore the patterns that 
occur when the terms of linear second-order recurrence sequences are 
combined in certain ways. First you will be asked to consider two different 
formulas involving the terms of the Fibonacci sequence, namely 


Pa 

—— and Fi+ Fi, 

and to look for patterns in the sequences defined by these formulas. You 
will then be asked to explore how these patterns generalise to other linear 


second-order recurrence sequences. 


Mathcad file 221A1-02 has been set up to help you to carry out this 
investigation — the worksheet displays the sequences defined by the above 
formulas. However, the notation u,, is used rather than F,, in the 
worksheet, because you can use the same worksheet to see the sequences 
corresponding to any linear second-order recurrence sequence u,, (subject 
to the accuracy of Mathcad calculations). You can do this just by 
changing the values of the initial terms a and }, and the coefficients p 
and q, in the recurrence system 


Uo = 4, U1 = b, Un+2 = PUn+1 a Gun (n ame 0, i, 2, me sds 


which defines the sequence u,, in the file. 


Mathcad notes provide extra 
information about the 
features and techniques used 
in the Mathcad files. They 


are optional. 


On the same line and to the 
right counts as ‘below’, but 
on the same line and to the 
left counts as ‘above’. 


Mathcad displays u2 as (un)?. 


The solutions use the 
notation F;,, since the 
sequence here zs the 
Fibonacci sequence, but you 


may choose to use u,, instead. 


If you try to display a table 
which contains a value that 
cannot be calculated, then 
Mathcad registers an error 
and shows no values at all. 


See Chapter Al, Section 3. 


CHAPTER Al 


Activity 4.3 Patterns in the Fibonacci sequence 


Open Mathcad file 221A1-02 Linear second-order recurrence 
sequences. The worksheet is set up with p= 1, q=1, a=0 and b=1, so 
u,, is the Fibonacci sequence. It displays the sequences 

Un+1 


and us +ur... 
Un 


n, Un; 
Look at each of the last two sequences in turn. For each sequence, try to 
spot a pattern, and hence make a conjecture about a general property of 
the Fibonacci sequence. 


Solutions are given on page 33. 


Mathcad notes 


© You may wonder why the range variable n used to display the tables is 
defined only to take values in the range 1,2,..., N — 1, when the 
recurrence system defines all the terms of the sequence from up to uy. 
Notice, however, that if we increase the final value of the range 
variable n from N — 1 to N, then the last term in the ‘ratio’ table 
should be uy4i/un. This term cannot be calculated, since the term 
Un+1 Was not defined earlier in the file. Also, if the first value of the 
range variable is 0, then the first term in the ‘ratio’ table should 
be u;/uo, which cannot be calculated, since Up = 0. 


The default result format has been set so that ‘Number of decimal 
places’ is 9. Also, ‘Exponential threshold’ has been set to 15, so 
numbers between 107!° and 10!° are shown in ordinary decimal 
notation. (Note that these settings affect only how Mathcad displays 
numbers — it always stores values internally to 15 significant figures for 
calculation purposes. ) 


© 


© All the tables have been resized, to show 19 values without scrolling 
(for the range 1,2,...,20— 1). If N is increased so that there are more 
values to display, then each table will scroll. 


In Activity 4.3, you saw that the sequence F,,,,/F;,, of ratios of successive 
terms of the Fibonacci sequence appears to tend to the golden ratio ¢, and 
to alternate above and below ¢. The golden ratio ¢ is one of the roots of 
the auxiliary equation associated with the Fibonacci sequence. 


In the next activity you are asked to use the same worksheet to explore the 
sequence Un+1/Un of ratios for a number of different linear second-order 
recurrence sequences u,,. You can change the recurrence sequence in the 
worksheet by altering the values of p, qg, a and 0. 


You will probably find it helpful to scroll down the worksheet until the line 
containing the definitions of p, g, a and b is at the top of your screen, 
before you begin Activity 4.4. This should allow you to see enough values 
of the sequence to spot any patterns, without having to scroll up and down 
every time you re-define p, q, a and 0. 


SECTION 4 _ Exploring linear second-order recurrence sequences with the computer 


Activity 4.4 Exploring un. /uUn 


For each set of values of p, g, a and 6 suggested in ‘Table 4.1, look at the 
sequence U,y+;/U, and make a note of any pattern you observe. Note that 
the file displays the roots of the auxiliary equation associated with the 
sequence u,, — which should be helpful! There are some spaces in the table 
for you to choose your own values of a and b, and you may wish to extend 
the table on a separate sheet to try further values of p and gq. 


Note that the table of ratios may register a ‘Found a singularity ...’ error 
for some combinations of p, g, a and b. This indicates an attempt to divide 
by zero, which happens here if one of the terms uj, Ug,...,Un_, of the 
original sequence is zero. 


When you have completed the table, try to make a conjecture about the 
long-term behaviour of the sequence Un+1/Un. 


Table 4.1 
Coefficients | Initial values Apparent pattern in 
p q| a b Unit han N= Lees 
1 1 ) 1 | tends to ¢, alternates above and below @ 
1 ae 1|tendsto® : 
1 1 
—1 an 1|teras to -@ 
= 1/1 1 teardD  - p 
—] 1 
3 reEra> Dd 
3 
3 


Solutions are given on page 33. 


In Activity 4.3 you met the conjecture that the sequence F? + F?,, 
consists of every other term of the Fibonacci sequence. In the next activity 
you are asked to explore the sequence u*, + uz,, for a number of different 
linear second-order recurrence sequences U,,. 


Activity 4.5 Exploring u2, + u2,, 
For each set of values of p, g, a and 6 suggested in ‘Table 4.2 overleaf, look 
at the sequence u? + u?,,, and make a note of any pattern you observe. 


When you have completed the table, try to make a conjecture about the 
sequence u? + u*,, which is more general than the above conjecture about 
Fe + F?,, in cases where p = 1, g = 1 and either a = 0 or b= 0. 


If you have time, you may also like to look at what happens if you vary p. 


You should still be working 
with Mathcad file 221A1-02. 
Ignore the sequence 

u2 + u2,, for the moment. 


If an error occurs, then 
Mathcad highlights the 
offending expression in red. 
Clicking on this expression 
reveals an error message. 


You should still be working 
with Mathcad file 221A1-02. 


CHAPTER Al 


Table 4.2 

Coefficients | Initial values Apparent pattern in 

p q\a b i te ag aed 2 Byi27) 

1 1/0 1 | it is every other term of uy, starting with u3 
l 1/0 2 

1 1 | 0 2 

1 ee: 0 

1 Live 0 

1 La 3 0 


Solutions are given on page 33. 


Now close Mathcad file 221A 1-02. 


4.3 Fibonacci numbers in sunflowers (Optional) 


Do not be tempted to spend The Mathcad file associated with this subsection allows you to investigate 
too long on this subsection! the close association between the pattern of florets on a sunflower head 
and the Fibonacci numbers. 


See Chapter A1, You may remember that the florets of a sunflower develop from tiny lumps 

Subsection 2.2. called primordia, which are created one after another on the edge of a disc 
in the centre of the sunflower head and then move outwards. In most 
sunflowers, the angle between successive primordia is about 137.5°. This is 
very close to 360(1 — 1/¢)°, where ¢@ is the golden ratio. 


The Mathcad file simulates the process that takes place on a sunflower 
head, and the pattern that it produces is shown in Figures 4.1 and 4.2. 
The file allows you to see what happens if you choose angles other than 
137.5°. It also allows you to investigate the spirals that appear in the 
sunflower pattern. Some clockwise spirals are highlighted in Figure 4.1, 
and some anticlockwise spirals are highlighted in Figure 4.2. We say here 
that a spiral is clockwise if your finger moves clockwise when tracing it 
inwards towards the centre, and anticlockwise otherwise. 


Figure 4.1 Clockwise spirals Figure 4.2 Anticlockwise spirals 
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SECTION 4 _ Exploring linear second-order recurrence sequences with the computer 


Activity 4.6 A Fibonacci sunflower (Optional) 


Open Mathcad file 221A1-03 A Fibonacci sunflower, and experiment 
in the ways described in the worksheet. 


Comment 


It seems that most choices of the angle between successive points produce 
a less well-packed pattern than that produced by the angle 137.5°. 


In the second part of the investigation, the angle between successive points 
is left fixed as 137.5°. It appears that highlighting every fth point 
produces a visible spiral if f is a Fibonacci number (other than 1, 2 or 3), 
but that no such spiral is produced for most other numbers f. 


It seems also that as f increases through the Fibonacci sequence, the 
spirals produced are alternately clockwise and anticlockwise, and 
progressively less ‘curly’. 


Now close Mathcad file 221A1-03. 


It is indeed true that Fibonacci numbers always give rise to spirals in a 
137.5° sunflower pattern in the way that we have described, and that these 
spirals have the properties mentioned in the comment above. If you are 
interested in the reasons behind these facts, you may like to read the 
explanation given at the end of this subsection. 


For each Fibonacci number f with f > 5, the spiral shown in the Mathcad 
file is just one of a set of spirals which do not intersect each other but 
together contain all the points in the sunflower pattern. The number of 
spirals in the set is f. In any particular sunflower pattern, just some of 
these sets of spirals are obvious to the eye — the ones in which successive 
points in each spiral are close together. 


You may be wondering how commonly real sunflowers display the type of 
pattern that we have explored. Research has established that 
approximately 95% of sunflowers are ‘Fibonacci’, where the angle between 
successive primordia is about 137.5°. Approximately 5% are ‘Lucas’, where 
the angle between successive primordia is about 99.5°, and the numbers of 
spirals are Lucas numbers. A small proportion of sunflowers do not fall 
into either of these two categories. Research also shows that successive 
primordia may appear in clockwise or anticlockwise order; 50% of 
sunflowers are clockwise and 50% anticlockwise. 


Explanation of the spirals observations 


The key to proving that the observations in the Comment following 
Activity 4.6 are true in general turns out to be the identity 


Pi = er ee = Oh 8. (4.2) 
which can be derived from Binet’s formula, as you will see in Section 5. 


The sunflower pattern in the Mathcad file is produced by plotting 500 
points. The first point is the one on the extreme right-hand side, and each 
subsequent point is plotted at an angle of approximately 360(1 — 1/¢) 
clockwise of the previous point, and a small distance nearer the centre. 


Figures 4.1 and 4.2 show some 
of the spirals corresponding 
te f = 2Zhandsf = 34. 


This botanical information 
was kindly provided by 

Prof. Ralph O. Erickson from 
the University of 
Pennsylvania. 


The Lucas sequence is 
2, 1,3, 4, 7,14, P2347, 76... 
See Chapter Al, Exercise 3.2. 


All angles in this discussion 
are measured in degrees. 


ui 


n=—I 
a Fy x 360 
ae tai eee 

: oe | —137.508 
ao, 2 84.984 
4 3 —52.523 
oS 32.461 
Ge -s —20.062 
(aoe 12.399 
o-cat —7.663 

9 34 4.736 
10 55 —2.927 


Remember that the first 
highlighted point is the one 
furthest from the centre of 
the sunflower pattern. 
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CHAPTER Al 


Since 


360 (1- 5) = 360 — ae 
«) w) 


the anticlockwise angle from each point to the next point in the sunflower 
pattern is 360/¢ ~ 222.5. Suppose that we choose the Fibonacci number 
F,,, and highlight every F,,th point. Then the anticlockwise angle from 
each highlighted point to the next highlighted point is 
360 

‘‘ 


By equation (4.2) (with n replaced by n — 1), this angle is equal to 
360 


Pr, % 


m1 
C) ae ee wy*) x =(F,...-% 360)-+ (= x 30) 
The first term in the expression on the right, F;,_; x 360, corresponds to 
F,,_; complete turns, and so the anticlockwise angle from each highlighted 
point to the next is just 


yr" 
p 


Now w = —0.618..., so w"—* alternates in sign and tends to 0 as n tends 
to infinity. It follows that the above sequence of angles behaves in the 
same way. Its terms for n = 2,3,...,10 are shown to three decimal places 
in the table in the margin. 


x 360. 


Thus, if n is even and large enough, then each highlighted point after the 
first is the same small number of degrees clockwise of the preceding 
highlighted point. Since it is also slightly closer to the centre, a clockwise 
spiral is produced. This is illustrated in Figure 4.1, which includes the 
spiral starting from the first point in the pattern when n = 8. There is a 
similar spiral starting from each of the first Fg points, making a set of 

Fs = 21 clockwise spirals in all, and three of these are shown in the figure. 
Similarly, if n is odd and large enough, then each highlighted point after 
the first is the same small number of degrees anticlockwise of the preceding 
highlighted point, and an anticlockwise spiral is produced. Figure 4.2, 
where n = 9, shows three of the set of Fy = 34 anticlockwise spirals. 


So for all large enough Fibonacci numbers F,, there are F,, spirals (provided 
that there are enough points in the sunflower pattern!), and the spirals are 
clockwise or anticlockwise according as n is even or odd. Also, since the 
magnitude of the angle between successive highlighted points decreases as 
n increases, bigger Fibonacci numbers produce less ‘curly’ spirals. 


Chapter A2, Section 6 
Conics on the computer 


In this section you will learn how to plot conics on the computer, using 
parametric representations. 


At the end of the section, there is an optional subsection in which you are 
invited to use Mathcad to explore how the shape of a conic changes as its 
eccentricity changes. 


6.1 Ellipses 


The standard parametrisation of the ellipse is See Chapter A2, 
P—ecost, y= being (8 < ¢ < 27). pene enernr 
This is a parametrisation of the ellipse in standard position with equation 
P+ ha (where a > b> 0). 
In the first activity you will plot ellipses in standard position, and a 
translated ellipse. 


Activity 6.1 Plotting ellipses from parametric equations 


Open Mathcad file 221A2-01 Ellipses in parametric form. Page 2 of 
the worksheet describes some basic Mathcad techniques for producing These techniques are used in 


graphs. the computer work for 


‘ : Chapter A2 of MST121, so 
When you have finished with page 2, work through pages 3 and 4 of the eemy be teoniliar with 


worksheet, carrying out Tasks 1 and 2. them already. 


Comments are provided in the worksheet, and there is a further comment 
below. 


Comment 


When you enter expressions involving a range variable in the z- and y-axis 
placeholders of a graph, Mathcad plots one point for each value of the 
range. By default, Mathcad joins these points together with line segments, 
to produce a continuous curve. The result is an approximation to the true 
curve — the more points plotted, the more accurate the approximation. 


As a rule of thumb, the range variable used to plot a curve should have a 
step size which ensures that at least 100 points are plotted. For the ellipses, 
the graph range t := 0,0.01..27 (a step size of 0.01) was used, which gives 
over 600 points. Note that the smaller you make the step size, the longer 
Mathcad will take to plot the graph, and once you reduce the step size 
beyond a certain level there will be no noticeable improvement in the 
shape of the curve. It is usually not worth plotting more than 1000 points. 
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Remember that Mathcad 
notes are optional. 


You should still be working 
with Mathcad file 221A2-01. 


See Chapter A2, 
Subsection 5.1. 


14 


CHAPTER A2 


Mathcad notes 


The technique used to plot two curves on a graph (such as the two ellipses 
on page 4 of the worksheet) can be extended to plot three, four, or more 
curves on the same graph. All the expressions for x should be entered in 
the x-axis placeholders, separated by commas, and the expresssions for y 
should be entered likewise in the y-axis placeholders. Mathcad plots the 
first y-axis expression against the first x-axis expression, the second 
against the second, and so on. When the graph is plotted, the line style 
and colour used to display each curve appear underneath the expressions 
on the y-axis. 


When a curve is represented by parametric equations, it is sometimes 
useful to be able to identify the point that corresponds to a particular 
value of the parameter. In the next activity you will see how to use 
Mathcad to identify a point on an ellipse. 


Activity 6.2 Identifying a particular point on an ellipse 


Move to page 5 of the worksheet and carry out Task 3. 


A solution and comments are provided on page 6 of the worksheet, and 
there are two further comments below. 


Comment 


© If you enter expressions involving an ordinary, single-valued variable 
(instead of a range variable) in the z- and y-placeholders of a graph, 
then the result is a plot of a single point. The point can be seen only 
after the graph trace has been formatted to display it as a symbol. 
This is the technique that was used to identify a particular point on 
the ellipse. 


© As the value of T increases through the range 0 < T < 27, the point 
(8cosT,5sinT’) travels once anticlockwise around the ellipse, starting 
and finishing at (8,0). In particular, T = 7/2 corresponds to the point 
(0,5), T = 7 to (—8,0), and T = 37/2 to (0, —5). 


Now close Mathcad file 221A2-01. 


6.2 Parabolas 


The standard parametrisation of the parabola is 

x is ae", ey ee Bat. 
This is a parametrisation of the parabola in standard position with 
equation 

y? =4ax (where a > 0). 


In Activity 6.3 you are asked to use Mathcad to plot such a parabola, with 
a= 1. There is no prepared Mathcad file for this activity — you are asked 
to create your own Mathcad worksheet. As with every worksheet that you 
create, it is a good idea to include some text, to make it more 
comprehensible to a reader such as your tutor (or yourself at a later date!). 


SECTION 6 Conics on the computer 


Activity 6.3 Plotting a parabola from parametric equations 


The instructions below tell you how to create a Mathcad worksheet 
containing a plot of the parabola with parametrisation x = t?, y = 2t. 


Part (b) describes how to enter a title in the worksheet, and you should 
also enter any other text that you think is appropriate. For example, when 
you define a variable, it is helpful to include some text nearby which 
indicates what the variable represents. 

(a) Begin by creating a new worksheet, as follows. Select the File menu 
and choose New.... In the list of templates that appears, Normal 
should be selected by default. If not, click on it. Then click on the 
OK button to create a new (Normal) worksheet. (Alternatively, type 
[Ctrl]n, or click on the New button on the standard toolbar.) 

(b) Enter a title at the top of your worksheet. To do this, first select the 
Insert menu and then choose Text Region. (Alternatively, type a 
double-quote ", given by [Shift]2.) Then type a suitable title — for 
example, Graph of parabola — in the text box. To finish, click 
anywhere outside the text box or press [Ctrl] [Shift] [Enter]. If you 
need to edit the text later, simply click on it. 

(c) Define a range variable t going from —5 to 5 in steps of 0.1. You can 
use the buttons on the ‘Calculator’ and ‘Matrix’ toolbars to do this, or 
just type t:-5,-4.9;5. Remember that the second number in the 
definition is the ‘next value’ in the range, not the step size. 

(d) Create your graph, positioning it below the definition of the range 
variable. To plot the parabola, enter t? in the x-axis placeholder (for 
example, type t A\2) and 2t in the y-axis placeholder (type 2*t). 


Fix the scales of your graph from 0 to 20 on the z-axis, and from —10 
to 10 on the y-axis, and resize the graph appropriately. 


(e) Select the File menu and use Save As... to name and save your 
worksheet. (You will be working on it again in Activity 6.4.) 


Comment 
A solution is shown below. 


Graph of parabola 


(sraph range t:=-5,-49.5 


10 


at 0 


-10 


Figure 6.1 The parabola y? = 4x in Mathcad 


See A Guide to Mathcad if 
you require more details on 
creating and editing your own 
worksheets. 


If you have just started 
Mathcad running, then there 
is no need to do this, as it 
automatically starts with a 
new (Normal) worksheet. 


Before you enter each new 
item in the worksheet, you 
should position the red cross 
cursor in an appropriate 
place. However, anything 
that you enter can be moved, 
or deleted, if you wish. 


You can click on the button 
on the ‘Graph’ toolbar, or 
type @, or use the Insert 
menu, Graph p X-Y Plot. 
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You should still be working 
with the Mathcad file that 
you created in Activity 6.3. 
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After the graph scales have been fixed as described in part (d) of the 
activity, the graph is 20 units wide and 20 units high. It should 
therefore be resized to make the graph box square, as shown above. 


When the axis scales are fixed as described, there are a few points on 
the parabola that Mathcad calculates but does not plot. For example, 
for the final value in the graph range, t = 5, Mathcad calculates the 
point (25,10). This point lies beyond the right-hand edge of the graph, 
so it is not plotted. When you use Mathcad to plot a graph, it does not 
matter if some points lie outside the graph box in this way. However, it 
is preferable (if not always possible!) to choose a combination of graph 
range and axis scales which avoids the calculation of a large number of 
points that lie outside the graph box. You should also try to avoid the 
calculation of points that lie a long way outside the graph box. 


In the next activity you are asked to identify a particular point on the 
parabola that you plotted in Activity 6.3. 


Activity 6.4 Identifying a particular point on a parabola 


Identify a particular point on the parabola x = i! y = 2tc-by plotting a 
second trace consisting of a single point (T*, 27’), as follows. 


(a) 


(d) 


First define a variable T with value —4; make sure that this definition 
is positioned above your graph. 


You may wish to make space for this new definition. To do so, position 
the red cross cursor where you want to insert some blank lines into the 
worksheet, then press the [Enter] key repeatedly until the space 
created is adequate for your needs. (Each key press will insert one 


blank line.) 


Position the vertical editing line at the right-hand end of the 
expression t? on the z-axis of your graph, with the horizontal editing 
line under the whole expression. (One way to do this is to click 
anywhere on the expression, then press [Space] and [Insert], as 
necessary, to select it all.) Type a comma, then enter T’ (for example, 
type ,TA2 — make sure that you type a capital T here). 


In the same way, position the vertical editing line at the right-hand 
end of the expression 2t on the y-axis and type a comma, then 
enter 2T (type ,2*T). 


Format the second trace to display the identified point as a box 
symbol. (Click in the graph to select it, and choose Graph » X-Y 
Plot... from the Format menu (or just double-click in the middle of 
the graph) to bring up the ‘Formatting Currently Selected X-Y Plot’ 
option box. Then choose the ‘Traces’ tab; then click on ‘trace 2’ and 
set the ‘Symbol’ to box and ‘Type’ to points. Click on the OK 
button to finish. ) 


Save your worksheet. 


You may now like to increase the value of T and check that the point 
(T*,2T) appears where you expect. 


SECTION 6 Conics on the computer 


Comment 


A solution is shown below. 
Graph of parabola 
Graph range t= -5,-49..5 


Particular value of the parametert T-=-4 
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Figure 6.2 Identifying a particular point on y? = 4x 


As T increases through the negative numbers towards 0, the point (T?, 2T) 
moves along the lower arm of the parabola towards the origin; when T = 0 
it is at the origin; and as T increases through the positive numbers it 
moves away from the origin along the upper arm of the parabola. 


If you found that you chose a value for T which caused the identified point 
to disappear from the graph, this was probably because the point lay 
outside the graph box! You can display its coordinates by evaluating T? 
and 27’ in the worksheet. You can see more of the parabola, and see the 
identified point for more values of J’, by changing the graph range and axis 
scales appropriately. 


Now close the Mathcad file that you have created. If you have made unsaved 
changes then you will see a 
dialogue box that asks 

6.3 Hyperbo las whether you wish to save the 
changes. 

The standard parametrisation of the hyperbola is 


ct=asect, y=btant (—ja1<t< ja,in<t< =a) 


This is a parametrisation of the hyperbola in standard position with 


equation 
2 y” 
i 1 (where a,b > 0). 


In the next activity you are asked to plot such a hyperbola, with a = 1 
and 6 = 2. Plotting a hyperbola is a little more tricky than plotting an 
ellipse or a parabola, because you have to plot the two branches as 
separate traces. 


You will use one range variable, tr, to plot the right-hand branch of the 
hyperbola, which corresponds to the range —}7 < t < }m, and a second 
range variable, t/, to plot the left-hand branch of the hyperbola, which 


corresponds to the range 1 + ST. 
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Activity 6.5 Plotting a hyperbola from parametric equations 


Open Mathcad file 221A2-02 Hyperbolas in parametric form, and 
carry out Task 1 on page 2 of the worksheet. 


Comments are provided in the worksheet, and there are further comments 
below. 


Comment 


© The graph ranges tr := —1.55, —1.54..1.55 and tl := 1.59, 1.60.. 4.69 
were chosen to ensure that no value of tr or tl is very close to —7/2, 
m/2 or 3/2. This avoids the calculation of points with very large 
coordinates. | 


© You may think that the initial graph shown in the comments in the 
Mathcad worksheet looks rather strange. This arises from the fact 
that if you do not enter numbers in the axis limit placeholders when 
you plot a graph, then by default Mathcad chooses limits which ensure 
that all the points that it calculates are plotted; it also chooses them 
to be appropriately round numbers. In this case, the endpoints of the 
portion of the hyperbola corresponding to the ranges of tr and tl are 


(sec(—1.55), 2 tan(—1.55)) ~ (48.1, —96.2), 
(sec 1.55, 2 tan 1.55) ~ (48.1, 96.2), 
(sec 1.59, 2tan 1.59) ~ (—52.1, —104.1), 
(sec 4.69, 2 tan 4.69) ~ (—44.7, 89.3). 
When the axis limits are set closer to zero, as suggested in the 


worksheet, the part of the curve near the origin can be seen more 
clearly, and the shape appears much more familiar! 


Mathcad notes 


The two-letter names tr and tl were used for the range variables because 
the ‘r’ and ‘I’ suggest their use in plotting the right and left branches of 
the hyperbola. However, we could equally well have used single-letter 
names, such as ¢ and u, or names consisting of a letter followed by a 
number, such as tl and t2. Note that care is needed when you use the 
letter ‘l’ in Mathcad, as it looks very like the number ‘1’! In general, 
Mathcad variable names can be any combination of letters and numbers, 
but they must start with a letter. 


aT 


In the next activity you are asked to add graph axes and asymptotes to 
the hyperbola that you plotted in Activity 6.5. 


The equations of the axes are z = 0 and y = 0. These can be plotted using 
the parametric equations 


g=0 yos end s=s-9= 6, 
respectively. 


The equations of the asymptotes of the hyperbola that you plotted in 
Activity 6.5 are y = 2x and y = —2z. These can be plotted using the 
parametric equations 


pee as an ee Se, ae, 


respectively. 


SECTION 6 Conics on the computer 


Activity 6.6 Adding graph axes and asymptotes 


Move to page 3 of the worksheet and carry out Task 2. You should still be working 


ith Mathcad file 221A2-02. 
Comments are provided in the worksheet. ~ ee ee 


Mathcad notes 


© There are four ways to add axes to a Mathcad graph: you can plot 
extra lines, format the graph to turn ‘Show Markers’ on, format the 
graph to turn ‘Grid Lines’ on, or format the graph to use the ‘Crossed’ 
axis style. You used the first of these in this activity. Details of the 
other three methods can be found in A Guide to Mathcad. 


© Note that the lines plotted for the axes may not align exactly with the 
tick marks for 0 at the edges of the graph box. The computer screen is 
made up of individual pixels (dots) which affect how accurately 
graphical information can be displayed, and there may be a tiny gap 
(of one or two pixels) between the lines and tick marks here. 


In the next activity you will explore how the shape of the hyperbola 
changes as you vary a and b. 


Activity 6.7 Exploring the shape of the hyperbola 


Move to page 4 of the Mathcad worksheet. This page contains a plot of You should still be working 
the hyperbola with parametrisation with Mathcad file 221A2-02. 


s=esect, y= Stent (—{4 <t< 50, o9 <t < Sm) 


together with its asymptotes, and the graph axes. It is set up so that you 
can vary the values of a and 6. 


(a) Keep a fixed at the value 1, and try the values 1, 2 and 3 for b. 
Describe the effect on the hyperbola. 


(b) Keep 6 fixed at the value 1, and try the values 1, 2 and 3 for a. 
Describe the effect on the hyperbola. 


(c) Explain the effects that you saw in parts (a) and (b). 


(d) Try the following pairs of values for a and b: a= 1, b= 2, and a = 2, 
b = 4. What is the relationship between the shapes of the two 
hyperbolas? You may find it helpful to print them out so that you can 
compare them. ; 


The point with parameter t = T is identified on the graph; you may like to 
try varying the value of TJ’ as well. 


Solutions are given on page 34. 


Comment 


As T increases through the range —}a < T < $n, the point 

(asec T, btanT) moves upwards along the right-hand branch of the 
hyperbola, lying on the x-axis when T = 0. Similarly, as T' increases 
through the range $7 < T < $7, the point (asecT,btanT) moves upwards 
along the left-hand branch of the hyperbola, lying on the z-axis when 

T =n. Values of T close to —}37, $m or $7 correspond to points that do 


~2 
not lie within the graph box. 
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Mathcad notes 


© You may have noticed that both graph scales are fixed from —S to S, 
where the variable S is defined at the top of page 4 of the Mathcad 
worksheet. This has the advantage that the graph scales can be 
changed just by changing the value of S — there is no need for any 
changes to the graph itself. Moreover, there is no need to change the 
range variable s which is used as the graph range to plot the axes and 
the asymptotes. This is because s is defined as s := —S,—-S+0.1..S 
(that is, s ranges from —S to S in steps of size 0.1), so the values 
taken by s will automatically be updated if S changes, and the axes 
and asymptotes will always be drawn to the the edges of the graph 
box. This is a useful general technique. 


© When a Mathcad graph is to include several curves, it is often best to 
plot them in the order ‘least important first’. This is done for the 
graph that you used in this activity — the axes are plotted as traces 1 
and 2, the asymptotes as traces 3 and 4, the branches of the hyperbola 
as traces 5 and 6, and the identified point as trace 7. The reason is 
that Mathcad draws trace 1 first, then trace 2 over the top of that, 
and so on. So a later trace may obscure some of an earlier trace. In 
fact, the trace order of the graph that you used in this activity does 
not have much effect, as the branches of the hyperbola cross the z-axis 
only once, and never cross the asymptotes! 


Now close Mathcad file 221A 2-02. 


6.4 Focus, directrix and eccentricity (Optional) 


The final, optional, activity in this section invites you to explore the shape 
of a conic with a fixed focus and directrix, and variable eccentricity. In 
particular, you can see the conic change between the three different types, 
ellipse, parabola and hyperbola, as the value of the eccentricity changes. 
The Mathcad worksheet also allows you to explore the focus—directrix 
property of the conics that it plots. 


Activity 6.8 Focus, directrix and eccentricity (Optional) 


Open Mathcad file 221A2-03 Focus directrix and eccentricity. The 
worksheet shows a plot of the conic with focus (1,0), directrix « = —1, and 
eccentricity e. The eccentricity e is initially set to 0.8, but you can vary its 
value. 


A point P on the conic is identified by a blue box symbol; you can change 
its position by changing the value of the variable T’. T'wo line segments 
meet at P; the length of the horizontal segment is the distance Pd of P 
from the directrix d, and the length of the other segment is the distance 
PF of P from the focus F’. The two distances PF and Pd, and the 

ratio PF'/Pd, are evaluated at the bottom right of the worksheet. 


(a) Vary the value of the eccentricity e, and observe how the shape of the 
conic changes. Some possible values for e are suggested in the 
worksheet, and you may like to try others of your own. In particular, 
you may like to explore the range of values of e for which it is difficult 
to distinguish by eye whether the conic is an ellipse, hyperbola or 
parabola. 


SECTION 6 Conics on the computer 


You may wish to look at a ‘small’ conic in more detail, or to see more 
of a conic that extends outside the graph box. You can decrease or 
increase the part of the plane corresponding to the graph box by 
changing the value of the variable S. 


Vary the value of T' (keeping e fixed) and observe how, although the 
distances PF and Pd vary as the point P moves along the conic, their 
ratio PF'/Pd remains constant — equal to the eccentricity e. 


Note that not every value of JT’ corresponds to a point on the conic. In 
fact, the x-coordinate of the identified point is equal to T’, so by 
looking at the graph scales you can see roughly which values of T’ 
correspond to points on the conic. If you choose a value of TJ’ that does 
not correspond to a point on the conic, then no point is plotted. 


Note that the graph is set up to identify points P on the upper half of 
the conic only. 


Comment 


When e is not much greater than 0, the conic is an ellipse, nearly circular 
in shape. It is positioned to the right of the y-axis, with two vertices on 
the x-axis. As e increases through the range 0 < e < 1, both of these 
vertices move along the x-axis: the left one moves further left and 
approaches the origin, while the right one moves further right, with the 
result that the ellipse becomes increasingly elongated. When e = 1, the 
vertex that was approaching the origin has reached it, while the other 
vertex has ‘disappeared to infinity’. ‘The conic is now a parabola in 
standard position. As e increases through the range e > 1, the vertex near 
the origin continues to move left along the x-axis, away from the origin, 
while the other vertex has ‘reappeared’ on the negative part of the x-axis, 
and moves right, approaching the first vertex. The conic is a hyperbola, 
whose asymptotes become increasingly steep. 


If you are interested in how the graph in the worksheet associated with 
Activity 6.8 is achieved, then you may like to read the following 
explanation. Otherwise, just close Mathcad file 221A2-03. 


Explanation of the Mathcad worksheet 


The worksheet states that the conic with focus (1,0), directrix x = —1 and 
eccentricity e is represented by the equation 


an” = 9° + 260 a0, (6.1) 


where a = e? —1 and B= e’* +1. 


If we then set x = t and solve equation (6.1) for y, we obtain 

y= +/a(t? + 1) +26t. Thus we can plot the conic on a Mathcad graph 
using two separate traces: the first with x-coordinate t and y-coordinate 
a(t? +1) + 26t, and the second with x-coordinate t and y-coordinate 
—/a(t? + 1) +26t. These are traces 5 and 6 in the Mathcad worksheet. 
Now t is defined to range from —S to S (where S is defined earlier in the 
worksheet), but not all of these values of t correspond to points on the 
conic. 


You may notice that a small 
gap appears in some cases 
between the upper and lower 
halves of the conic. This is 
explained, along with other 
aspects of the worksheet, at 
the end of this subsection. 


For example, with e = 0.8, 
you can see that values of T 
between about 0 and 9 
correspond to points on the 
conic (the exact range 

is ¢ <T <9). 


This equation can be 
obtained using the method of 
the solution to Chapter A2, 
Activity 3.5. 


a1 


Note that 


Sa faBy. : 
<7 


This assumes that S = 10 in 
the worksheet, so that t takes 
the values 


10. 9.98.4 40. 


Le 
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So what happens for the other values of t? In fact, t corresponds to a point 
on the conic if and only if a(t? + 1) + 2G¢ is non-negative. If a(t? + 1) + 2(¢ 
is negative, then according to its definition, the y-coordinate of the point 

to be plotted is the square root of a negative number. This does not make 
sense, and so Mathcad does not plot a corresponding point on the graph. ° 


The coordinates of the identified point P on the graph are x = T, 

y = \/a(T? + 1) + 267, so the identified point always lies on trace 5, the 
upper trace. If T is set to a value that does not correspond to a point on 
the conic, then no point is plotted. 


We end with explanations of three other aspects of the graph. 


© The line segments indicating the distances PF’ and Pd on the graph 
are drawn by plotting the three points (uo, vo) = (1,0), (ui, v1) = (2, y) 
and (u2,v2) = (—1,y), with a line trace to join them together. This is 
trace 7 on the graph. The calculations are off the Mathcad page to the 
right. 


© The graph range is defined as t := —S,—0.998S ..S; that is, t ranges 
from —S to S in steps of 0.0025. This gives a way of plotting 1000 
points, whatever the value chosen for S. 


© In some cases a point where the conic intersects the x-axis does not 
correspond with a value of t in the chosen graph range. Then Mathcad 
cannot plot the corresponding point on the x-axis, which may lead to 
a slight gap between the upper and lower traces for the conic close to 
that point. 


For example, with e = 2.7, the hyperbola cuts the z-axis at (—#, }s 


but —i — —0.459... is not a value taken by t within the graph range. 
The closest such values are t = —0.46, for which no point is plotted, 
and t = —0.44, for which the corresponding points on the conic are 
approximately (—0.44, +0.46). Hence a gap of height 2 x 0.46 = 0.92 


appears on the Mathcad graph. 


Now close Mathcad file 291A2-08. 


Chapter A3, Section 5 
Isometries on the computer 


In this section, you will see how Mathcad can demonstrate the effect of 
applying isometries in R*. First isometries are applied to triangles. Then 
Mathcad is used to plot the graphs of conics with ry-terms, by applying 
isometries to simpler conics. 


There is also an optional subsection in which you are invited to use 


Mathcad to explore surface and contour plots of functions of two variables. 


5.1 Isometries and triangles 


In this subsection you will use Mathcad to explore the effect of isometries 
on triangles. The notation used in the Mathcad worksheet for this 
subsection differs in some respects from that used in the earlier sections of 
Chapter A3. These changes are needed to allow isometries to be 
implemented and composed effectively in Mathcad. 


The main difference is that points in R° are represented using vector 
notation: for our current purposes, a vector is a column of numbers. For 
example, to define the point P with coordinates (5,—2) as a vector in 
Mathcad, we use 


px (5). 


You can think of P as representing a sequence of two numbers Py) = 5 and 
P, = —2, written in a vertical table. Thus P is the x-coordinate of P and 
P, is the y-coordinate of P. 


(5.1) 


We define an isometry in Mathcad as a function whose inputs and outputs 
are vectors. For example, to define the translation t that moves each point 
two units to the right and one unit down, we use 


ial Cen 


Once this definition has been made in a Mathcad worksheet, the function t 
can be applied to any input vector. For example, if P and t have been 


defined as in equations (5.1) and (5.2), then evaluating t(P) gives a 


(5.2) 


Also, evaluating t(P)o gives 7 and evaluating t(P), gives —3. 


Vector notation is used throughout the Mathcad worksheet for this 
subsection, but you will not need to edit any vectors in the worksheet, nor 
to create any new ones. 


In the first activity you will see vector notation used to implement a 
translation in Mathcad and demonstrate its effect on a triangle. 


A vector is a particular type 
of matrix (a rectangular 
array of numbers). You will 
study matrices in detail in 
Block B. 


Remember that in Mathcad 
the first term of a sequence 
has subscript 0, by default. 


In Chapter A3, Section 2, we 
used t2,-; for this translation, 
but in Mathcad this notation 
is not convenient. 


You will learn how to enter 
and edit vectors, and matrices 
in general, in Block B. 
Instructions are also given in 


A Guide to Mathcad. 
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Similar advice applies to 
many of the Mathcad graphs 
that you will see in this 
section. See the solution to 
Activity 5.5, for example. 


Remember that when plotting 
the effects of isometries it is 
important to have equal 
scales on both axes. 
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Activity 5.1 Points, triangles and the translation function 


Open Mathcad file 221A3-01 Isometries and triangles, and move to 
page 2 of the worksheet. Read this page carefully — it introduces the 
notation that will be used throughout the worksheet. 


To check that you understand the way in which the translation function 
has been implemented, try changing the values of c and d so that the 
image triangle is in a different position, for example, with its left-most 
vertex at the origin. 


Comment 


© The purpose of the ‘point function’ 


Pin (i.) 


defined in the worksheet, is to give an easy way to deal with the 
vertices (0, Yo), (£1, Y1), (2, y2) and (x3, y3) = (Xo, Yo) of the triangle. 
It means that P(n) is the nth vertex, and t(P(n)) is the image of the 
nth vertex under the translation t. The x-coordinate of the nth vertex 
is then P(n)o, and its y-coordinate is P(n);. Similarly, the 
z-coordinate of the image of the nth vertex under t is t(P(n))o, and its 
y-coordinate is t(P(n));. You can see these expressions used on the 
axes of the graph to plot the two triangles. 


© Inthe Mathcad graph, the original triangle is plotted as a solid black 
trace and the image triangle as solid red. If you have difficulty 
distinguishing between the triangles or wish to print the graph ona 
non-colour printer, then you may prefer to change the line style of the 
image triangle trace from ‘solid’ to ‘dot’ or ‘dash’. 


Mathcad notes 


The graph at the bottom of the Mathcad page has been resized and both 
scales have been fixed from 0 to 10. It has also been formatted to show 
erid lines — this was done by switching on ‘Grid Lines’, then switching off 
‘Auto Grid’ and setting the ‘Number of Grids’ to 10 for each axis. 


The Mathcad notation that you saw in Activity 5.1 is used throughout the 
rest of the Mathcad worksheet to implement isometries and demonstrate 
their effects on triangles. 


In the next activity you will see the effect of a rotation demonstrated in 
Mathcad. The rotation is implemented by defining an appropriate value 
for the angle @ and setting 
(Py Po cos(6) — P, sin(@) 
MAT = \ P. sin(@) + P, cos(9) }’ 


Remember that P, and P, are the x- and y-coordinates of P. 


SECTION 5  Isometries on the computer 


Activity 5.2 Rotations 


Move to page 3 of the Mathcad worksheet. A triangle is defined near the 
top of the page. Further down, the rotation r = rg is defined, and a graph 
shows the effect of rg on the triangle when 0 = ST. You can change the 
value of 6, and you can enter this angle in either radians or degrees, as 
explained in the worksheet. 


(a) Figure 5.1 shows the original triangle defined in the worksheet, and its 
images under three rotations, namely r3, jay Tonja O00 T_5,76.. Try to 
decide which image corresponds to which rotation. Then confirm your 
answer by trying each of the three values of 6 in turn in the worksheet 
and checking that the image appears where you expect. 


Figure 5.1 


You should still be working 
with Mathcad file 221A3-01. 


2) 
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(b) Repeat part (a) for Figure 5.2, which shows the images of the triangle 


under rotations through 60°, —160° and —100*°. 


Figure 5.2 


Solutions are given on page 34. 


Comment 


Note that traces 1 and 2 on the graph plot the z- and y-axes, using the 
parametrisations r= s, y=Oandx=0,y=s. 


Mathcad notes 


© 


Note that deg (in lower-case letters) is a built-in Mathcad constant. 
Its value is 0.01745..., the number of radians in one degree. For 
example, evaluating 180 deg gives 3.14159.... (The constant deg is 
included in the system of units built into Mathcad, alongside standard 
SI units, such as m for metre and kg for kilogram. So deg 1s not 
defined in worksheets where these units are switched off - Math 
menu, Options... , ‘Unit System’ tab, Default Units ‘None’.) 


You may find it helpful to reformat the graph to add grid lines. To do 
this, switch ‘Grid Lines’ on, switch ‘Auto Grid’ off, and set ‘Number of 
Grids’ to 20 for both axes. (You can also choose the ‘Grid Color...’ — 
the default colour for the grid lines is green.) However, displaying this 
number of grid lines makes the axis labelling rather overcrowded, and 
may make the graph less clear if printed. 


o_O 


In the next activity you will see the effect of a reflection demonstrated in 


Mathcad. 


SECTION 5 Isometries on the computer 


Activity 5.3 Reflections 


Move to page 4 of the Mathcad worksheet. This page is set out ina similar You should still be working 


way to page 3. The reflection q = qg is defined, and a graph shows the with Mathcad file 221A3-01. 

effect of gq, on the specified triangle, together with the line of reflection. We use @ for the angle of a 

You can change the value of ¢. rotation and ¢ for the angle 

(a) Figure 5.3 shows the original triangle defined in the worksheet, and its oes to avoid 
confusion. 


images under three reflections, namely qs5x/6, d3x/4 and q3x jg. Lry to 
decide which image corresponds to which reflection. Then confirm 
your answer by trying each of the three values of ¢ in turn in the 
worksheet and checking that the image appears where you expect. 


Figure 5.3 


(b) Edit the Mathcad page to change the original triangle to the triangle 
with vertices (0.9,6.2), (5.7,2.6) and (6.3, 8.4), and set ¢ to 45°. Then, Recall that you can enter an 
by experimenting with different values of ¢, find to the nearest degree angle in degrees by using the 


the angle ¢ such that the corresponding reflection maps the triangle built-in Mathcad constant 
onto itself. ‘deg’. For example, for an 


angle of 180°, enter 180 deg 


Solutions are given on page 34. (type 180*deg). 

Comment 

The line of reflection is plotted using the parametrisation The line making the angle ¢ 
with the positive x-axis has 


f=. Be stand. slope tan ¢. 


In the next activity you will see the effect of composites of isometries 
demonstrated in Mathcad. 


27 


You should still be working 
with Mathcad file 221A3-01. 
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Activity 5.4 Composite isometries 


Move to page 5 of the Mathcad worksheet. This page demonstrates the 
effect of the composite go f, of an isometry f followed by an isometry g, 
on a triangle. Initially, f is defined to be the reflection qy where ¢ = 0, and 
g to be the translation t..g where c = 1 andd=0,sogof isa 
glide-reflection parallel to the z-axis. You can see different glide-reflections 
by changing the values of ¢, c and d appropriately (where ¢ is the angle 
that the line y = (d/c)x makes with the positive a-axis). 


Figure 5.4 shows the original triangle defined in the worksheet, and its 
images under three glide-reflections, namely to,2 ° dz/2, t-6,-6 © In/4 

and t_11,.0 ° go. Try to decide which image corresponds to which 
elide-reflection. Then confirm your answer by trying each of the three sets 
of values of ¢, c and d in turn in the worksheet and checking that the final 
image appears where you expect. 


Figure 5.4 


Solutions are given on page 34. 


Comment 


If you wish, you can experiment with the composites of different types of 
isometries, by changing the definitions of f and g. 


i 


Now close Mathcad file 221A3-01. 


SECTION 5 Isometries on the computer 


5.2 Isometries and conics 


In the main text, we developed a method of sketching a quadratic curve L 
with an zy-term. This method involves first determining a conic K with 
no ry-term that is congruent to L, and a rotation that maps K onto L, 
then drawing a graph of K, and finally applying the rotation to obtain a 
graph of L. The final two stages of this process can be carried out using 
Mathcad, as you will see in this subsection. 


The Mathcad worksheet for this subsection uses vector notation in a 
similar way to the worksheet for Subsection 5.1. 


In the first activity you are asked to plot a quadratic curve that is the 
image under a rotation of an ellipse in reflected standard position. You 
will use a page of a Mathcad worksheet that has been set up to plot the 
ellipse with equation 

2 2 


i Y 
a ae (a,b > 0), (5.3) 
and its image under a rotation rg, for values of a, b and @ that are defined 
in the worksheet and that you can change. 


Activity 5.5 A rotated ellipse 


Open Mathcad file 221A3-02 Isometries and conics, read the 
introduction on page 1 of the worksheet, and move to page 2. 


On this page, the ellipse with equation (5.3), and three isometries t = t..4, 
r=frg and q = qj, are defined. An isometry f is initially set to be the 
rotation r, and the original ellipse and its image under f are plotted on a 
graph. 


You can change the values of a and 6 to specify different ellipses. You can 
also change the definition of f to see the image of the ellipse under 
different types of isometries, and you can change the values of c, d, 9 and ¢ 
to change the isometries themselves. 


The quadratic curve L with equation 
192? + 6ry + 1ly? — 40 =0 
is the image of the ellipse, in reflected standard position, with equation 
5u* + ty? = 1 
under a rotation of approximately 18° about the origin. Change the values 
of a, 6 and 6 on page 2 of the worksheet, to obtain the graph of L. 


A solution is given on page 34. 


Comment 


The original ellipse is plotted using the technique that you saw in the 
computer work for Chapter A2: a large number of points on it are 
calculated using a suitable range variable u and the usual parametrisation 
of the ellipse (which remains valid for an ellipse in reflected standard 
position). However, here we use vector notation to make it easy to deal 


See Chapter A3, 
Subsection 4.2. 


This ellipse is in standard 
position if a > 6, and in 
reflected standard position 
if b> a. 


See Chapter A3, Example 4.1. 


Remember that rotations 
through a positive angle are 
anticlockwise. 


In Chapter A2 the name t 
was used for the range 
variable, but here we use wu. 
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You should still be working 
with Mathcad file 221A3-02. 


You could work this out for 
yourself by following the 
strategy given in Chapter A3, 
Subsection 4.2. 


30 


CHAPTER A3 


with the points on the ellipse and their images under the isometry J. oust 
as the ‘point function’ P defined the vertices of the triangle in 
Subsection 5.1, so here we have a ‘parametrisation function’ p, given by 


__ (acos(u) 
4d Pam ( retire 
which defines points on the ellipse. Each value of u corresponds to a point 
p(u) on the ellipse, and the image of this point under the isometry ‘i 
is f(p(u)). You can see these expressions used on the graph axes to plot 


the ellipse and its image under f; as usual, the subscripts 0 and 1 specify 
z- and y-coordinates, respectively. 


Mathcad notes 


© The square root operator can be obtained by clicking on the 
appropriate button on the ‘Calculator’ toolbar, or by typing \ 
(a backslash). For example, to obtain /2 you can type \2. 


© You cannot use the same name for a variable and a function in a 
Mathcad worksheet. Here we use t¢ for a translation, so a different 
name is required for the range variable — we use uw. 


In the next activity you will use Mathcad to produce a graph of a 
hyperbola whose equation has an ry-term. 


nnn nEEEEEEE Ea 


Activity 5.6 A rotated hyperbola 


Move to page 3 of the Mathcad worksheet. It is set up in the same way as 
page 2, but the curve plotted is a hyperbola. 


The quadratic curve L with equation 
24x? + 20xry — 24y” — 13 =0 

is the image of the hyperbola with equation 
24" — Ge? =4 


under a rotation of approximately 11° about the origin. Change the values 
of a, b and @ on page 3 of the worksheet, to obtain the graph of L. 


A solution is given on page 34. 


a 


We now discuss quadratic curves that are the image under a rotation of a 
hyperbola in reflected standard position. 


A hyperbola in reflected standard position does not have parametric 
equations of the the same form as those for a hyperbola in standard 
position, so the original hyperbola on the Mathcad page will never be in 
reflected standard position. 


Instead, we can start with a hyperbola in standard position, reflect it in 
the line y = x to obtain a hyperbola in reflected standard position, and 
then rotate it through an appropriate angle to obtain a graph of the 
required quadratic curve. You are asked to do this in the next activity. 


SECTION 5 Isometries on the computer 


Activity 5.7 A reflected and rotated hyperbola 


The quadratic curve L with equation 
x” + 12ry + 6y? — 30=0 


is the image of the hyperbola K in reflected standard position with 
equation 
ge 2 
2 ge We 
Mo <8 
under a rotation of approximately —34° about the origin. This hyperbola is 
in turn the image of the hyperbola H in standard position with equation 
fe 2 
ee, 
ae 


under reflection in the line y = z, that is, under qz/4. 


(a) On page 3 of the Mathcad worksheet, set 6 = 0 (if it is not already set 
to 0), and change the values of a and b to obtain a graph of H. The 
original and image hyperbolas will then be superimposed. 


In the definition of the isometry f, change r to q; this means that f is 
the isometry gy. Then change the value of ¢ to obtain a graph of the 
hyperbola K. 


(c) Now edit the definition of the isometry f again so it reads as follows: 
f(P) = r(q(P)). 


This means that f is the isometry rg o gy. Change the value of 0 to 
obtain a graph of L. 


(b) 


Solutions are given on page 35. 


Comment 


Since reflecting a hyperbola in standard position in the line y = x produces 
the same image as rotating it about the origin through the angle 90°, an 
alternative way to obtain a graph of the above quadratic curve L is to 
begin with the same hyperbola H in standard position, and simply rotate 
it through 90° + (—34)° = 56°. 


In the next activity you will use Mathcad to produce a graph of a parabola 
with an zy-term. Here, we need to use a translation as well as a rotation. 


Activity 5.8 A translated and rotated parabola 


Move to page 4 of the Mathcad worksheet. It is set up in the same way as 
pages 2 and 3, but the curve plotted is a parabola and the isometry f is 
set initially as a translation. 


The quadratic curve L with equation 
9r* — 24ry + 16y? — 10x — 70y — 75 = 0, 

is the image of the quadratic curve K with equation 
y’ — 2x — 2y —3 =0, 


which has no zy-term, under a rotation of approximately 37° about the 
origin. Completing the square in the equation for K, and simplifying, gives 


(y — 1) — 2(n + 2) =0, 


See Chapter A3, Activity 4.3. 


Reflecting in the line y = x 
corresponds to exchanging 
the roles of the x- and y-axes. 


You should still be working 
with Mathcad file 221A3-02. 


Note that in order to obtain 
the graph of L, you only need 
to set the values of a, 6, ¢ and 
0, and edit the definition of f 
as in part (c). 


You should still be working 
with Mathcad file 221A3-02. 


You could work this out for 
yourself by following the 
method explained in 
Chapter A3, Subsection 4.2. 
Note that in this case D 4 0 
and FE £0. 
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See Chapter A3, Section 1.3. 


This is the ellipse EF’ discussed 
in Chapter A3, Section 1. 


32 


CHAPTER A3 


so the curve K is itself the image of the parabola H in standard position 
with equation 


y = 2a, 
under the translation t_1. 


(a) On page 4 of the worksheet, set c = 0 and d = 0 (if they are not 
already set to 0) and change the value of a to obtain a graph of H. 


(b) Change the values of c and d to obtain a graph of K. 
(c) Now edit the definition of the isometry f so that it reads as follows: 
f(P) = r(t(P)). 


This means that f is the isometry rg 0 t..4. Change the value of @ to 
obtain the graph of L. 


Solutions are given on page 35. 


Now close Mathcad file 221A8-02. 


5.3 Surface and contour plots (Optional) 


In this subsection you are invited to use Mathcad to explore surface and 
contour plots of functions of two variables. A contour plot can provide a 
useful check on the shape and position of a quadratic curve. 


Activity 5.9 Surface and contour plots (Optional) 


Open Mathcad file 221A3-03 Surface and contour plots, and read 
through pages 1 and 2 of the worksheet. 


(a) In the definition of the function f on page 2 of the worksheet, change 
‘d’ to ‘g’ so that the function g is plotted instead of the distance 
function d. 


The quadratic curve with equation 
An” + 9y? — 162 — 18y —11=0 
is then displayed as the contour with height 0. 


(b) Now use page 2 of the worksheet to verify the shapes of the quadratic 
curves that you considered in Activities 5.5 to 5.8, as follows. 


For each quadratic curve, change the values of the coefficients in the 
worksheet to the values of A, B, C, D, E and F corresponding to the 
equation of the curve. Check that the shape and position of the 
contour with height 0 are the same as those of the curve. 


(c) If you would like to see how to create and to format surface and 
contour plots, then read page 3 of the worksheet. 


a 


Now close Mathcad file 221A3-03. 


Solutions to Activities 
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Solution 4.3 


The value of the expression F;,,,/F, appears to tend 
to the golden ratio @ as n increases, and appears to 
lie alternately above and below @. 


(You met this property in Section 2 and will see a 
proof of it in Section 5.) 


The value of the expression F? + F?,, always seems 
to be a Fibonacci number. These numbers appear in 
alternate positions in the F,, table. To be precise, 


F24+F2 = Fy, 
F3 + F2 — Fs, 
2 = 


This suggests the conjecture that, for n = 1,2,3,..., 
| ake oe Jae = Fon41- 


(This conjecture is true. It can be proved using 
Binet’s formula, but we shall not do so.) 


Solution 4.4 


It appears that the sequence u,+41/un always tends 
to one of the roots of the auxiliary equation, provided 
that these roots are real. In fact, if the roots are real 
and distinct, then the sequence seems to tend to the 
root of larger magnitude, whether this is positive or 
negative. (The magnitude of a real number z is z if 
x > 0 and —z if x < 0. For example, the magnitudes 
of 3 and —3 are both 3.) It is natural to conjecture 
that this happens for any linear second-order 
recurrence sequence whose auxiliary equation has 
real roots. Changing the initial values a and b does 
not seem to affect this long-term behaviour. 


(This conjecture is indeed true for all such linear 
second-order recurrence sequences u,,, but we prove 
it only for the particular case of the Fibonacci 
sequence, in Section 5.) 


You may also have noticed that although for the 
Fibonacci sequence the terms of the sequence 
Un+41/Un appear to alternate above and below the 
value to which they tend, this does not happen for 
every linear second-order recurrence sequence. 


(This behaviour is related to whether the roots of the 
auxiliary equation are of the same sign or of opposite 
signs. ) 


Solution 4.5 


It appears that if p= 1, q= 1 and a = 0, then the 
sequence u2 + u2,, (n = 1,2,3,...) consists of b 
times every other term of the original sequence u,, 
starting with us. This suggests the conjecture that if 
p=1,¢q=1 and a =, then, for alfd, 


ae Pe bee pis for’ n= 15 2) 3; (1) 


Similarly, it appears that if p= 1, q=1 and d= 0, 
then the sequence u% + u2,, (n = 1,2,3,...) consists 
of a times every other term of the original sequence, 
starting with u2. This suggests the conjecture that if 


p=1;¢=1 and 6= 0, then, for all a, 
i, see (2) 


If you tried varying p, then you would have found 
that this seems to have no effect on the above 
conjectures; that is, the first conjecture seems to be 
true if g= 1, a= 0 and for all p and b, and the 
second seems to be true if g = 1, b= 0 and for all p 
and a. 


(The identities in equations (1) and (2) can in fact be 


deduced from the identity 
as + Fe, Sag OP = 12.8, «.., 


which you met in Activity 4.3. For equation (1), we 
use the fact that if p= 1, q = 1 and a = 0, then 


20). et aS 1s 
4 beretore, for 7 =. b,2.3, ...« , 
Un + Ungar = (OF)? + (bFn41)? 
= 6° (Fe + Fi, 4) 
Se Pens i 
== Dilan a1. 


Equation (2) can be proved in a similar way, using 
thetact. that afp = 1, q¢= 1 and 6= 0, then 


e255 for mn =, 2,3,.... 


In fact, there is a general identity 
qu + ers Spree legs 4, for ti =:¥, 2, 3,....4 


which holds for all values of p, g, a and b. This can 
be proved using the closed form of the sequence, but 
the details are a little involved.) 


= fe 


COMPUTER BOOK A 


Chapter A2 


Solution 6.7 


(a) Ifa is fixed, then as b increases, the two points 
where the hyperbola crosses the x-axis remain 
fixed, and the asymptotes become steeper. 


(b) If b is fixed, then as a increases, the two points 
where the hyperbola crosses the z-axis move 
away from the origin, and the asymptotes 
become less steep. 


(c) The points where the hyperbola crosses the 
x-axis are (ta,0), and the asymptotes are 
y = +(b/a)x, where a,b > 0 (see Section 2). 
Thus if a is fixed, then the two crossing points 
remain fixed; if b now increases, then b/a 
increases and hence the asymptotes become 
steeper. On the other hand, if b is fixed and a 
increases, then the crossing points move away 


from the origin; also b/a decreases and hence the 


asymptotes become less steep. 


(d) The two hyperbolas have the same asymptotes, 
since they have the same value of b/a. In fact, 
their shapes are even more closely related: they 
are similar. The first hyperbola has 
parametrisation 


=> sect, y= ane 


3 
(—in<t< 5m,50<t< 57) 


and the second has parametrisation 


x's Zeect,. y= 4tant 


1 3 
(—in<t< $n,50<t< 5n), 


so if the point (u,v) lies on the first hyperbola, 
then the point (2u, 2v) lies on the second. 


(In general, any two hyperbolas with the same 
value of b/a are similar: each is obtained from 


the other by a scaling with the same scale factor 


with respect to both axes. This is also true for 
ellipses. ) 
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Solution 5.2 


(a) A is the image of the triangle under Pan j4, Bis 
its image under r_s,/¢6, and C is its image 
under T_/4- 


(b) A is the image of the triangle under a rotation 
through 60°, B is its image under a rotation 
through —160°, and C is its image under a 
rotation through —100°. 


34 


Solution 5.3 


(a) A is the image of the triangle under q3,/g, B is 
its image under q3,/4, and C is its image 
under q5x/6- 

(b) Reflection in the line through the origin that 
makes an angle of approximately 53° with the 
positive x-axis maps the triangle onto itself. 


Solution 5.4 


A is the image of the triangle under to,2 © g,/2, B is 
its image under t_¢,—6 © dx/4, and C is its image 
under t_j1,0 © qo- 


Solution 5.5 


Setting a = 2, b = 2 and 6 = 18° gives the graph in 
Figure $3.1, in which L is the dashed trace. Here the 
axes limits have been changed to make the plot of 
the ellipses larger. 


3 
2 
| Ps 
: ! 
“a 
“S “2 a | 0 1 2 


“3 


Lad 


Figure S3.1 


Solution 5.6 


Setting a = 1/./2, b= 1/2 and 6 = 11° gives the 
graph in Figure $3.2, in which L is the dashed trace. 


NZ 
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Figure S3.2 


Solution 5.7 


(a) Setting 6=0, a = V3 and b = V10 gives the 
graph of H, as in Figure $3.3. 


4 
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-4 “3 0 2 4 
Figure S3.3 


(b) Setting ¢ = 45° gives the graph in Figure S3.4, 
in which K is the dashed trace. 


mee 
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Figure S3.4 


(c) Setting 06 = —34° gives the graph in Figure $3.5, 


in which L is the dashed trace. 
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SOLUTIONS TO ACTIVITIES 


Solution 5.8 


(a) Setting c= 0,d=0 and a= § gives the graph 
of H, as in Figure 83.6. 
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Figure S3.6 


(b) Setting c= —2 and d= 1 gives the graph in 
Figure $3.7, in which K is the dashed trace. 


-4 = 0 2 4 
Figure $3.7 


(c) Setting 0 = 37° gives the graph in Figure S3.8, 
in which L is the dashed trace. 


-4 4 0 2 4 
Figure S3.8 
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